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Abstract 

We present a general method for the derivation of various statistical 
quantities describing the detection of a beam of atoms emerging from a 
micromaser. The use of non-normalized conditioned density operators 
and a linear master equation for the dynamics between detection events 
is discussed as are the counting statistics, sequence statistics, and waiting 
time statistics. In particular, we derive expressions for the mean number 
of successive detections of atoms in one of any two orthogonal states of 
the two-level atom. We also derive expressions for the mean waiting times 
between detections. We show that the mean waiting times between de- 
tections of atoms in like states are equivalent to the mean waiting times 
calculated from the uncorrelated steady state detection rates, though like 
atoms are indeed correlated. The mean waiting times between detections 
of atoms in unlike states exhibit correlations. We evaluate the expressions 
for various detector efficiencies using numerical integration, reporting re- 
sults for the standard micromaser arrangement in which the cavity is 
pumped by excited atoms and the excitation levels of the emerging atoms 
are measured. In addition, the atomic inversion and the Fano-Mandel 
function for the detection of de-excited atoms is calculated for compari- 
son to the experimental results of Weidinger et al. jl| , which reports the 
first observation of trapping states. 



1 Introduction 



The most fundamental system for studying matter-radiation coupling is a single 
two-level atom interacting with a single mode of an electromagnetic field — first 
treated theoretically by Jaynes and Cummings Q . The one-atom- maser or mi- 
cromaser Q is an experimental realization of this system. Other applications 
for the Jaynes-Cummings model include ion traps Q, quantum nondcmolition 
measurements ^ 0] , quantum state teleportation H, ^ , quantum computa- 
tion jl^ , and optical cavity QED . For a nice review of the early experiments, 
see Raithel et al. jl^. For more recent discussions, consult [l5[ |. 

Much work has been done on determining just how much information about 
the micromaser cavity field can be obtained from the detection statistics of 
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the emerging atoms. In the standard micromaser experiments, atoms enter the 
high-Q resonator in the upper of the two Rydberg states of the maser transition 
and the emerging atoms are probed for being in either of these two states. For 
this case, the detection statistics of the emerging atoms are sensitive only to the 
photon-number distribution in steady state. The detection statistics have been 
shown to contain information on the photon statistics and the intensity 
fluctuations 0, ^ |§. In order to measure field coherence properties, 
means for breaking the phase symmetry of the field have been proposed — both 
by coherently pumping the cavit y l2l|,M p3| and by a measurement induced 
breaking of the phase symmetry p.4\, 25|, 26|. Both arrangements can be used 
to measure the phase diffusion rate and lincwidth of the field, but neither has 
been shown to give the full details of the spectrum. Cresser has recently 
shown that all of the information needed to obtain the full spectrum is contained 
in the atomic beam and has proposed a Mach-Zehnder type interferometer for 
measuring the general micromaser spectrum. 

Most of the efforts to make connections between the detection statistics and 
the field properties has concentrated on the correlation functions for the de- 
tection of de-excited atoms, which are not conditioned to exclude detections 
between events. Knowledge of the whole set of correlation functions would pro- 
vide a complete description of the statistical properties from which all statistical 
quantities of interest could be calculated in principle; however, many statistical 
quantities (e.g. waiting times) are much easier to derive from the conditioned 
evolution of the cavity field between detections. 

The conditioned evolution of the cavity field between detections must in- 
clude the effect of the passage of atoms undetected due to detector inefficiency. 
As shown by Briegel et al. this leads to a nonlinear master equation for 
the conditioned evolution of the field between detections. However, a linear 
master equation and a non-normalized conditioned density operator, as used by 
Herzog [jl^ , can be used in place of the nonlinear equation and the normalized 
conditioned density operator. In Sec. |^ of this paper, we will show how the 
two methods are equivalent and review the methods of calculating detection 
statistics in the micromaser in general. 

Using a perturbation method based on the nonlinear equation of Ref. |^ , 
McGowan and Schicve have derived an effective micromaser spectrum as calcu- 
lated from the known detection statistics — one that agrees with the true spec- 
trum in the case of perfect detector efficiency, but otherwise includes a back- 
action effect due to measurement. Herzog points out that this is not the 
conventional micromaser spectrum; usually defined as the Fourier transform of 
the steady state normalized two-time correlation function of the electric field 
strength. The differences between the two definitions for the spectrum is an 
interesting question for the future that will not be discussed in this paper. 

In Sec. |, we discuss three types of statistical quantities: counting statis- 
tics, sequence statistics, and waiting time statistics. In Sec. 4.1 we will discuss 
counting statistics and, in particular, the Fano-Mandel function for measuring 
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the difference in the variance of the counts from that of Poissonian statistics. 
They were treated rather nicely by Briegel et al. ||2^ and will be discussed here 
only briefly. 

In Sec. 4.2, we will discuss sequence statistics. Sequence statistics are fun- 
damentally characterized by the probabilities for a given number of successive 
detection events being of atoms in a particular sequence of the states |A) and 
|B), where I A) and |B) are any two orthogonal states of the two- level atom. 
In Sec. |4.3| , they will be used to derive expressions for the mean number of 
successive detections of atoms in state |A) and the mean number of successive 
detections of atoms in state |B). The method used is analogous to the one used 
by Englert et al. to derive an expression for the mean number of successive 
detections of atoms in the same state. 



Waiting time statistics will be discussed in Sec. 4.4. There we will derive 



expressions for the mean waiting times between successive detections. We show 
that the mean waiting times between detections of atoms in the same state 
are equivalent to those calculated from the uncorrelated steady state detection 
rates. Of course atoms in the same state are correlated, and we show that the 
simplification that leads to the uncorrelated mean waiting times between atoms 
in the same state does not occur for mean squared waiting times and presumably 
for higher powers of waiting times between atoms in the same state. We also 
derive expressions for the mean waiting times between atoms in unlike states, 
which do exhibit correlations. 

In Sec. 4.5 , we will discuss a numerical integration method based on the linear 
master equation which can be used to evaluate the expressions derived in Sec. ^, 
and in Sec. |^ we report results for the standard micromaser arrangement. In 
addition to evaluating the expressions derived in Sec. |4| we calculate the atomic 
inversion for comparison to the new experimental results of Weidinger et al. ^ 
who recently reported the first observation of trapping states. We find that 
there is qualitative agreement. In an attempt to improve the agreement, we 
included an uncertainty in the interaction time with a Gaussian distribution, 
among other effects, but failed to improve the agreement significantly. We also 
evaluate the Fano-Mandel function for the detection of de-excited atoms for 
comparison to the results of Weidinger et al. We find poor agreement due in 
part to the level of amplitude suppression exhibited by the experimental results, 
presumably due to some effect not included in our calculation. We evaluate the 
mean number of successive detections of atoms in each state and the mean 
waiting times between detections of atoms in unlike states, and find that they 
are strongly affected by the presence of trapping states — offering an alternative 
means of the experimental observation of trapping states. 
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2 Micromaser Dynamics 



The density operator of the micromaser cavity field evolves as 



Pit) 



Xrt 



P(0) 



(1) 



where r is the rate of atomic injection with assumed Poissonian statistics, and 
the evolution operator X is given by 



X = C + M{ip) - 1 



(2) 



Here C is the Liouvillian operator that governs thermal damping of the field, 
and M.{f) is a super operator that governs the change in the field due to the 
passage of a single two-level atom with accumulated Rabi angle ip. 
The field damping operator L is given by [M 



Cp- 



l 



[{i' + 1) (a^ ap — 2apa) + pa^a) + v [aa^ p — 2a^ pa + paa^)] , (3) 



where N^yi = r/7 for convenience, v is the mean number of thermal photons, 
and 7 is the photon decay rate. Note that for any p, tr{£p} = 0. 

The atomic passage operator A^((p) results from solving the atom- field in- 
teraction problem for the combined atom-field density operator Pa-f and then 
tracing over the atomic subsystem. By performing the trace in a basis spanned 
by two orthogonal states |A) and |B) of the two-level atom, M. can be written 
S.S M = A + B where Ap = (A|pa-f |A) and Bp = (B|pa-f |B). The states |A) 
and |B) may be the Rydberg states of the maser transition, or any coherent 
superposition of those states. In this paper we will consider only the case in 
which all of the atoms arrive in the upper maser level. Then is given by [p^ 



Mp 



cos 



^(p\/aaJ^ p cos ^i^Voflt^ -f 7^ sin aat^ p sin ^(^\/aaJ^ 7 , (4) 



with the normalized raising and lowering operators 

I I 

7' 




(5) 



The form of M is different for the general case (in which the atoms arrive in a 
coherent superposition). Note that for any p, trjAtp} = tr{p}. 

The steady state density operator is defined by Xp^^ = 0. When all of the 
atoms arrive in the upper maser level {M given by Eq. (|4|)), p^^ is given by 
Ref. ISlll 



p-(ata)=p-(0) n 



(6) 
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where p^^{0) is determined by the normaUzation condition trjp''^} = 1. 

The accumulated Rabi angle is given hy ip — gtmt , where tint is the interaction- 
time and 

9^ J J dxg{x) , (7) 

is the atom-field coupling constant given in terms of the position dependent 
atom-field coupling strength g(x) for a cavity of length L. 

As described by Meystre et al. in a cavity field with no photons, each 

atom will undergo an integer number q of Rabi cycles if obeys the equation 

'^(""''^^^ ^ ^ 

When this equation is satisfied, each atom leaves the cavity photon number 
unchanged, and the photon number is "trapped" with uq being the maximum 
number of cavity photons allowed (except for thermal noise). These trapping 
states will be shown to strongly affect the statistical quantities derived in Sec. ^. 



3 Detection Statistics 

For a particular detection event of interest, it is convenient to define the operator 
corresponding to the occurrence of that detection event, and the operator 
= X — X^ corresponding to evolution of the cavity field in the absence of 

detection events of that type. For example, if we are interested in the detection 

of outgoing atoms in state |A), then X+ — ijaA and X~ — C+{1 — i]a)A+B — 1, 

where rjA denotes the detector efficiency. 

A detection event at time t leads to the non-local state reduction 

X+p{t) 

P^'^ ^ ir{x^p{t)} ' 

and occurs with a probability density 

rtY{X+ p{t)} . (10) 

For times between detections, the passing of an atom at time i, undetected 
due to detector inefficiency, leads to the non-local state reduction 

{M - X+)p{t) 
P^'^ tr{{M-X^)p{t)} ' 

and occurs with a probability density 

rtr{{M- X+)p{t)} . (12) 
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This, along with the continuous effect of thermal damping, leads to the master 
equation for the evolution of the density operator between detections condi- 
tioned by the absence of detections 

= A-pe-trlA-pelPc , (13) 

where we have used the identities ti{A4pc} ~ 1, = X — X~ , and, for any p, 
ti{Xp} = 0. This was derived by Briegel et al. @ for the case X^ = ryA-^+TyeS, 
corresponding to the detection of atoms in either state (Eq. (2.25) of that paper). 
This equation is nonlinear, and has the formal solution given by Eq. (2.26) of 
Ref. 11 

The nonlinearity in Eq. ( |l3|) is what leads to the normalization of the conditioned 
density operator (|T^). 

The conditional probability density for a detection occurring at time t given 
that no detections have occurred between times and t is given by 

rlY{X+pS)} ■ (15) 

The exclusion probability, or the probability that a detection does not occur 
between times and t is 



exp 
If we define 



ii(i)^tr{e'^"'-*p(0)} , (17) 

We find that 

u{t) = rtr{X-e^~''* p{Q)} ^ -rlviX+e^'"-* p{Q)} , (18) 

where we have used X^ = X — X^ and, for any p, trjA'p} 0. So the 
expression ( p^ becomes 

"■'■(/'"|j)^«<"' f^") 

Thus the exclusion probability is u{t) given by Eq. (^) and is equivalent to the 
denominator of Eq. dlj), as was noted by Englert et al. pO[. 
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The probability density for the next detection occurring at time t is the 
conditional probability density for a detection occurring at time t multiplied 
by the probability u{t) that the condition is satisfied, or 

u{t)rtr{X+p,{t)} - rtr{A'+pe(i)} . (20) 

Here we have defined the non-normalized conditioned density operator 

PeW^«WPcW=e^"'-*p(0) , (21) 

where the decay in the norm of the operator is due solely to the exclusion of 
detections between times and t. This is the solution to the linear master 
equation 

I'Pc = X'Pc ■ (22) 

This has been used by Herzog |Q . The two approaches to Eq. ( pO| ) are identical, 
but the linear equation (p^), as opposed to the nonlinear equation (p^, is simpler 
to numerically integrate (see Sec. 4.5 for a discussion of om numerical integration 
method). 

In the following sections, the subscripts "A", "B", and "AB" will denote 
detections of atoms in state |A), state |B), and state |A) or |B) (e.g. = 
ijaA + tibB). 



4 Observables of Interest 

Here we will discuss three types of statistical quantities: counting statistics, 
sequence statistics, and waiting time statistics. Counting statistics and the 
Fano-Mandel function were treated rather nicely by Briegel et al. and will be 
discussed in Sec 



and. In Sec. 4.3 



4.1 only briefly. Sequence statistics will be discussed in Sec. 4.2 



will be used to derive expressions for the mean number of 



successive detections. In Sec. 4.4, we will derive the mean waiting times between 
successive detections. In Sec. 4.5, we will discuss a numerical integration method 
which can be used to evaluate the expressions. 

Ignoring prior observations, detections of atoms in state |A) and |B) will 
occur at the uncorrelated rates 



(23) 



For any statistical variable x, it is interesting to calculate the average normalized 
to the uncorrelated value 



(24) 
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where (x) is the average of x, and (a;)uncor is the average performed in the 
absence of correlations with the detection rates (p3|). This ratio is equal to 
unity for Poissonian statistics, and otherwise tells us whether the correlations 
are positive or negative, and, for the case of successive detections, whether the 
detection events are bunched or antibunched. 

As the detectors become less efhcient, the times between detections grow 
larger (with an increasing number of atoms passing undetected between detec- 
tions). The detected atoms thus become less correlated and we have the general 
result 

for any x , (x)iiorin 1 , as 7?a, ??b ^ . (25) 
4.1 Counting statistics 

The counting statistics were discussed by Briegel et al. |28| . By counting statis- 
tics, we mean the statical properties of the detection events counted for a certain 
observation period. They are fundamentally characterized by the joint proba- 
bilities for detecting n atoms in state |A) and m atoms in state |B) in an ob- 
servation time t. Of particular interest is the Fano-Mandel function Q which 
measures the deviation of the variance of the counts from that for Poissonian 
statistics Q . It is defined by convention as 

Q{t) ^ ((7V(i)-(7V(i)))2)_-l, 
{N{tf) 



(Nit)) 



{N{t)) - 1 , (26) 



where N{t) is the number of counts in an observation time t. An approximate 
expression for calculating the Fano-Mandel function was developed by Rempe 
and Walther |]l6j , who achieved results that agreed well with a computer simula- 
tion, and later with experimental results | |3^ . An exact expression was derived 
by Briegel p8[ | , and was shown to agree well with the previous approximation in 
the regime for which it is valid. For counting atoms in state |B) in an observation 
time t, the Fano-Mandel function can be expressed as |^ 

Qb(0= 2r^ di' (^1-0 tr{A'+e-^'^*'XBP^^} , (27) 

where 

MbP^-^^-P- (28) 

The Fano-Mandel hmction Qa for the detection of atoms in state |A) can be 
obtained by replacing with X^ in Eqs. ( p7| ) and (|2^). Qb is related to 
the two-time correlation function for the detection of atoms in state IB) lESl. 
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The Fano-Mandel function QfioW for the number of photons in the cavity can 
be related to the short-time hmit of the same atomic correlation function Jist , 
providing a connection between Qb and Qfioid- 



4.2 Sequence statistics 

Here we ask a different question: "given n successive detections, what is the 
probability that the detections are of atoms in a particular sequence of the 
states |A) and |B)?" We use the notation P[ABBA. . . {n terms)] to denote the 
probability of an n-detection event being comprised of atoms detected in the 
sequence of states "ABBA . . . ." They are normalized such that, for each n, the 
sum of all P[x], where x is a sequence of "A" and "B" of length n, is unity. 
The simplest of these probabilities are the probabilities that a single atom is 
detected in state |A) or |B), given by 

At the next level of complexity are the probabilities P[AA], P[BB], P[AB], 
and P[BA] for n = 2. Rather than calculate them directly, it is convenient to 
first consider the conditional probabihties P[AA], P[BB], P[AB], and P[BA], 
where the underline indicates the given condition. For example, P[AB] denotes 
the probability that the second detection of a two-detection event is of an atom 
in state |B), given that the first detection is of an atom in state |A). The joint 
probability P[AB], for example, is then given by the conditional probability 
P[AB] multiplied by the probability P[A] that the condition is satisfied: 

P[AB] = P[A] P[AB] . (30) 

Note that this relationship between joint probability and conditional probability 
is true regardless of what condition is made. For example, it is also true that 
P[AB] = P[B]P[AB]. 

We will now derive an expression for P[AA]. Expressions for the other 
conditional probabilities can be generated in an analogous manner. Given that 
an atom in state |A) was detected at time we have p(0) Xj^ p'^ytriXj^ p^^}. 
Until a later time t, when a second atom in state |A) is detected, we whish 
to exclude detections of atoms in either state, so we use the non-normalized 
conditioned density operator pc{t) — e'^AB''*p(0). Then the probability density 
for the time t when the second atom is detected is rtr{X^pc{t)}, so that 

P[AA] = r r dti,{X+p,{t)} = ^ tr [x+-^X+pA . (31) 



The identities 



Xp'"" = , 

for any p , iY{X p} = , (32) 
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are especially useful. Using = X — Xj^q — X^ and the identities (^2|), we 
have 

tr lx+-±X+A = tr{X+} - r (33) 



and 



where we have defined 



rstrJ<-^A-+p'4 =tr|A-+-i<p''i. , (35) 



'■A y,- ^^B f ^ "M '^'■B ^^A I 

^AB J ^AB 

The sequence probability P[AA] and the other sequence probabilities for 
n = 2 are given by 

P[AB] = P[BA] = I . (38) 

Eq. ( |3^ ) provides us with an interpretation of F as relating to the probability 
that a switch occurs in the type of detection event, given by 

P[AB]+P[BA]^ . (39) 

H^abP } 

Note also that these probabilities obey a useful "distributive property". For 
example, 

P[AA] + P[AB] ^ P[A(A + B)] P[A] , (40) 

which is a way of stating that the detection of an atom in state |A) will nec- 
essarily be followed by the detection of an atom in either state |A) or |B) with 
unit probability. This distributive property holds for all sequence probabilities. 

Expressions for sequence probabilities for n = 3 and greater can be formed 
by a sequence of the operators X^ and X^ in right to left order, separated by 
the operator — l/A'^g. For example, 

1 , r,.^ -1 -1 



PfABB] = —-^ -tr <^ X+ X+ A-^p"" )■ , (41) 



stated without proof. 
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4.3 Successive detections 



The mean number of successive detections of atoms in the same state, (n), 



was first calculated by Wagner et al. |24, E5l E6| using Monte Carlo techniques 



for the phase-sensitive micromaser experiment described in those papers. Later, 
Englert et al. derived an expression for (n) and (n)norm- Using an eigenvalue 
method, they reproduced the results of Wagner et al. with good agreement and 
calculated (n.)norm for the standard micromaser experiment. Here we will derive 
the mean number (tt-a) of successive detections of atoms in state |A) and the 
mean number (ub) of successive detections of atoms in state |B). The derivation 
of (n) = |((riA) + ("-b)) in Ref. |Q needs to be modified only slightly in order 
to obtain (tia) and (tt-b). 

Given that the detection of an atom in state |B) has occurred, we denote 
the probability of detecting n atoms in state |A) prior to the next detection of 
an atom in state |B) by = P[BA"B], where "A"" is short for n terms of "A" 
in the sequence, (similarly for p^). They obey the normalization 



E 



oo 

P^ = J2Pn=^- (42) 

=0 



Note that the probability of having n detections of atoms in one state between 
two detections of atoms in the other state, p„, used in the derivation of (n) in 
Ref. is related to and pj^ by 

Pn = P[B] + P[A] p^ = P[BA"B] + P[AB"A] . (43) 

We then consider the probability = P[BB] of detecting an atom in state 
|B) after a previous atom detected in state |B), and the possible ways in which 
the two events can occur between the detection of two atoms in state |A): 

Po = P[BB] = P[ABBA] 

+ (P[ABBBA]+P[ABBBA]) + --- . (44) 

Noting that, for example, 

P[ABBA] ^ P[A] P[ABBA] = P[B] P[ABBA] , (45) 
Eq. (U) becomes 

P^ = ^(P[ABBA]+2P[ABBBA] + ...) 

^ ' n=2 

A similar expression is obtained for . Together with Eqs. (^), we have 

i-P^-^f^-P^ + ^-^y (47) 
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and a similar expression for p^. 

We now consider the probability P,^ of detecting n atoms in state |A) in 
succession (similarly for P^). They are normalized as 

oo oc 

T.Pn=T.Pn=l- (48) 

n—1 n—1 

differs from only by the exclusion of n = 0, so that 

For the mean number of successive detections of atoms in state |A) wc obtain 
K) ^ 5]nF„A = -^ = -tr{A'+p-}, (50) 

n—1 ^ 

and similarly 

("B> - - ^trjA-B+p-} . (51) 



The average of the two expressions gives 

1 1 



tr{A'+Bp-} , (52) 



' ' 2P[AB] 2r 

which is the inverse of the switch probability (^9|) and is equivalent to Eq. (3.24) 
of Ref . |P0| . The uncorrelated values are 



/"OO 

(p^)uncor - re / dt e-('-*+'--)* = P[B] 

JO 



(53) 

and similarly for (pq )uncor and (nB)uncor- The average of (nA)uncor and (nB)uncor 
gives Eq. (3.10) of Ref. |§ 

^ (54) 



^.^unco, 2P[A]P[B] ■ 

Thus 

_ P[A] P[B 



(,flA;norm = [nB 



" ^ P[AB] 

tr{.y+p--} tr{X+p^^} 
tr{A'+Bp-}r 



(55) 
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4.4 Waiting time statistics 

Waiting time distributions have been calculated previously for particular trap- 
ping states |2^, |l^. Here we will calculate the mean waiting times between 
various detection events. 

We begin by deriving an expression for the mean time (tA-+A) between de- 
tections of atoms in state |A). Given that an atom in state |A) was detected 
at time we have p{0) = X^^ p^^ /tr{X^ p^^} . At a later time t, the next atom 
in state |A) is detected. Until then, we whish to exclude detections of atoms in 
state I A), but we do not care how many atoms are detected in state |B). So we 
use the non-normalized conditioned density operator pc{t) — e'*A''*p(0). Then 
the probability density for the time t when the next atom is detected in state 
I A) is rti{X^pc{t)}. The mean time between detections of atoms in state |A) 
is then 



{tA^A) = r [ dtttr{X+p,{t)} = — tr J X+^X+p^"" \ 
Jo rA [ J 

Using X^ = X — X^ and the identities ( |3^ ) we have 



(56) 



trl^X+^X+p^^^=l, (57) 

and we conclude that the mean time between detections of atoms in state |A) 
is equal to the uncorrelated value 

f°° 1 
(iA^A> uncor = rA dtte-'^^* = —. (58) 

Jo ^A 

Likewise, for atoms in state |B), (tB^s) = (^B^B)uncor = l/^^B- 

These results do not prove that there are no correlations between atoms in 
the same state — only that we have to go to higher powers in the waiting time 
in order to see those correlations. For example, the mean squared time between 
detections of atoms in state |A) is 

+2 







(4^A> = r I dtt'tr{X+p,it)} 



rrA 1 




-H 

rrA 1 







(59) 

where, in the last step, we have again substituted X^ = X — X^ and used the 
identities (^2|). This does not reduce to the uncorrelated value 

2 2 

(*A-*a) uncor — T • (60) 

TA^ 
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If we want a more interesting mean waiting time (one that exhibits correla- 
tions), we have to compute (^a^b) or {tB~>A)- Starting with the same initial 
condition p(0) = A'/p''7tr{A'/p'^"}, and using the non-normalized conditioned 
density operator pc{t) — e'^B '"*p(0) for the absence of detections of atoms in 
state |B), the probability density for the time t until the next detection of an 
atom in state |B) is rtr{X^pc{t)}. The mean time until the next detection of 
an atom in state IB) is 



/■°° 1 f 1 
(U^b) = r / dtttr{X+pS)} ^—tv{ X+ :^X+p^^ 

Jo fA [ 

Using = X — X^ and the identities ( |3^ ) , we have 



so that 



and similarly for (^b^a)- The uncorrelated value is 



so that 



(61) 



tr \ X+^X+p^^ = tr{-iA'+p-} , (62) 



(tA^B>--tr{-iA'+p-}, (63) 



(iA^B)uncor = ^B / dite"'''^* = — , (64) 



(tA^B)norm = — tr{-^X + p''} , (65) 
TA X^ 



and similarly for (tB^A)- 

4.5 Numerical integration method 

Traces of the form 

^ = tr|0i^03P^^| (66) 
can be evaluated by first numerically integrating the linear equation ( [2^ ) 

'p,{t) = -02Pcit) , Pc(0)=O3P^^ , (67) 

to obtain Pc{i), and then integrating the equation 

X{t) - tr{Oipe(i)} , ^(0) - , (68) 
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to obtain X — linit^oo X{t). 

When the eigenvalues and eigenstates of the operator O2 are known ana- 
lytically, the eigenvalue method used in Ref. |^0| is much faster. However, for 
the general case, the eigenvalues and eigenstates must be obtained numerically. 
If one does not care about the eigenvalues, then numerical integration is sim- 
pler to implement and likely to be faster. In calculating the mean number of 
successive detections of atoms in the same state, the Authors of Ref. ^ state 
that numerical integration is faster "roughly by a factor of two" , but neither 
the numerical integration method nor the method for obtaining the eigenvalues 
and eigenstates were clearly stated. 

As an alternative to numerical methods, McGowan and Schieve have 
used a perturbation method to derive approximate analytical solutions for (n). 
This is a general method that could be used, in principle, to calculate other 
statistical quantities, though it is unclear as to when the method is a valid 
approximation. 

5 Standard Micromaser Experiment 

Excluding added complications (e.g. two atom events or non-Poissonian injec- 
tion statistics), the methods for deriving various statistical quantities describing 
the detection of the emerging atoms (Sec. ^ are completely general, as are the 
expressions given in Sec. For a particular micromaser experiment only the 
super operators A and B, and the steady state density operator p^^ need to be 
specified. Here we will give results only for the standard micromaser arrange- 
ment in which the atoms arrive in the upper maser level, and the excitation 
level of the emerging atoms are measured. Other micromaser experiments may 
be considered in future papers. 

We have allready specified p^^ (Eq. (^)). All that remains is to specify A 
and B. Identifying the states |A) and |B) with the upper and lower maser levels, 
respectively, A and B are given by Q 

Ap = cos{ipV aa'f)p cos{(fV aa^) , 

Bp = 7^ sin((/3Vaat)psin((y9\/aaJ)7 . (69) 

Note that only the diagonal elements of the operators A, B, and £ are needed, 
since p^^ is diagonal, and these operators preserve diagonality. 

All numerical results presented in this paper assumed symmetric detector 
efficiencies r?A = ??b = though the methods can be used for non-symmetric 
detector efficiencies. 

5.1 Atomic inversion 

Trapping states were recently observed by Weidinger et al. Q in measurements 
of the atomic inversion. In that paper, the atomic inversion is defined as the 
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difference in the probabilities of a given injected atom being detected in the 
ground state and the excited state 

/EEtr{(A-+-A-+)p-}-i(rB-rA) . (70) 

Another definition is the difference in the probabihties of a given detected atom 
being in the ground state and the excited state 

/ = P[B] - P[A] ^ . (71) 

Both formulas reduce to the difference in the number of atoms in the ground 
state and excited state tr{(S — A)p^^} in the ideal limit of one hundred percent 
detector efficiencies. 

Fig. 1^ shows / calculated for A'cx = 7 and Nc^ = 10 for comparison to 
Figs. 3(a) and 3(/3) of Ref. §. Here v = 0.054 and rj = 40%. The coupling 
constant for the 63P3/2^61I?5/2 Rydberg transition in rubidium was determined 
in Ref. |l] as 5 = 39 ± 5 kHz. Here we have used g = 39 kHz to plot / as 
a function of tint. The vertical dotted lines indicate the same trapping states 
as in Fig. 3 of Ref. [Q . It should be noted that the linear trend present in the 
raw experimental data due to the unequal spontaneous emission rates and its 
subsequent removal makes the scale on the vertical axis on Figs. 3(a) and 3{(3) 
of Ref. [0 useful only for comparisons of relative heights. There is qualitative 
agreement between the theoretical and experimental results in that each set 
of plots has dips centered on the trapping states and approximately the same 
amplitude of modulation, but little quantitative agreement. As a simple attempt 
to improve the agreement, we added the effect of an uncertainty to tint with a 
gaussian distribution to simulate velocity averaging in the experiment. The 
dashed curves on Fig. |] are the effect of an uncertainty of 3 fis. This has 
the effect of producing a level of amplitude suppression comparable to that for 
the experimental data. Other effects that we added (but do not include the 
results here) were the spontaneous emission of the atoms and the removal of 
the subsequent linear trend, and a misalignment of the detection unit (which 
would cause incorrect detection) as discussed in Sec. 3.1.4 of Ref. ||l^. The 
inclusion of either of these effects did not cause any significant change in the 
shape of the graphs. A more involved parameterized fitting of the experimental 
results correcting for the effects of detuning, spatial dependence of the atom- 
field coupling, and others (see the discussion of Eq. (73) in Ref. |jl^) may prove 
useful, but will not be attempted here. It should be emphasized that this 
lack of quantitative agreement is within the margins of statistical error of the 
measurements, and does not negate the observation of trapping states or indicate 
any problems with the theory. 
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5.2 Fano-Mandel functions 

Figs. H and ^ show the steady state Fano-Mandel functions Qa(oo) and Qb(oo), 
scaled by I/77A and l/jye respectively. The results are independent of detec- 
tor efficiencies when scaled in this way. The parameters used here and in all 
subsequent plots are N^^ = 7 and i/ = 0.054. Fig. || shows many regions of 
sub-Poissonian fluctuations in the detection of de-excited atoms, i' — 0.054 is 
enough to destroy the visibility of most trapping states (not shown) . For smaller 
values of trapping states have the effect of reducing the fluctuations of counts, 
sometimes to sub-Poissonian levels, creating dips in the Fano-Mandel functions. 
Exceptions occur for Qb(oo) at the trapping states (p = n, 27r, Stt, . . . ,, where 
maxima occur in the fluctuations. 

Fig. ^ shows Quit) for observation times of one and four cavity decay times, 
as well as an average of Qb {t) for observation times between one and four cavity 
decay times. This is to be compared to Fig. 4 of Ref . |^] , where the points plotted 
were averaged over time periods between one and four cavity decay times. Here, 
we again have rj = 40% and have used g = 39 kHz to plot Qb (t) as a function of 
tint- The vertical dotted lines indicate the same trapping states as in Fig. 4 of 
Ref. [0. In sub-Poissonian regions, the Fano-Mandel function reaches a steady 
state value at about one cavity decay time |2^, which can be seen here by 
the small amount of differences between QB{t = I/7) and Qb(^ = 4/7) in the 
sub-Poissonian regions. In the super-Poissonian regions, QB{t) approaches its 
steady state value at a much slower rate. 

Compared to the theoretical calculation, the experimental results exhibit an 
amplitude suppression of roughly a factor of ten. Sources of amplitude suppres- 
sion include velocity averaging and the spontaneous decay of the atoms during 
the times between leaving the cavity and reaching the detector. A discussion 
of these added effects and the resulting fluctuations may be included in a later 
paper, but will not be included here. Without the inclusion of these effects, 
the comparison may not be a meaningful indicator of the amount of agreement 
between the experimental results and the theory. 

5.3 Successive detections 

Fig. 1^ shows the mean number of successive detections of excited atoms, {ua), 
for T] — 100%, 40%, and 10%. The vertical dotted lines indicate the trapping 
states for no = 0, . . . , 4. At a trapping state, if it were not for the effect of 
thermal damping, we would expect all of the atoms to emerge in the upper 
maser level, having undergone an integer multiple of Rabi cycles. Some trapping 
states give rise to large peaks indicating an increase in the number of successive 
detections of excited atoms as expected. Others have a smaller effect on the 
shape of the graph that can only be seen in the change of the slope of a nearby 
peak. This suggests that the curve could be approximated by a superposition 
of suitably weighted gaussians centered on the trapping states. Fig. shows the 
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mean number of successive detections of de-excited atoms, (ne). This graph 
shows the analogous effect of trapping states with the trapping states giving 
rise to dips rather than peaks. Both Figs. || and ^ show a decrease in the 
mean number of successive detections as the detectors become less efficient. The 
uncorrelated values are unaffected by the detector efficiencies (for equal detector 
efficiencies), and this decrease in the mean number of successive detections is just 
an indication of the decrease in correlations for decreasing detector efficiencies. 
Fig. shows the mean number of successive detections (of any type) normalized 
to the uncorrelated value, (nA)norm = ("■B)norm = (''T')norm- For decreasing 
detector efficiency, (n) approaches unity as expected. Antibunching occurs in a 
narrow interval of = 9. 

5.4 Waiting times 

Fig. 1^ shows the average time until the next detection of a de-excited atom 
after an initial detection of an excited atom, (tA->B)j scaled by rrjB- The scaling 
compensates for the increase in mean waiting time as the detectors become elss 
efficient and the probability of the final detection occurring decreases. Again, 
the vertical dotted lines indicate the trapping states for uq = 0,...,4. The 
trapping states have the effect of giving rise to peaks in the average waiting 
time for the next detection of a de-excited atom as expected. As before, some 
of the trapping states have a small effect. Fig. ^ shows rrjAitB^A) which shows 
the analogous effect of trapping states giving rise to dips rather than peaks. 
Fig. and |l| show (tA^B)norm and (tB^A)norm respectively. Fig. |l^ shows 
some negative correlations in a limited vicinity of </? = 9. 



6 Conclusions 

The effect of trapping states on the mean number of successive detections of 
atoms in a particular state and the mean waiting times between detections of 
atoms in unlike states offer an alternative means of the experimental observa- 
tion of trapping states, though a comparison of the errors in calculating these 
quantities from experimental data is required to determine which is most useful. 
In the near future, we plan to consider the phase-sensitive setup of pof , 
in which case only the form of the operators A and B need to be altered for 
the expressions derived in this paper to carry over. We may consider other 
experimental setups as well. 

It would appear that the theory of atomic detection statistics in the micro- 
maser is mature, though it would be interesting to include effects absent in the 
treatment of the ideal micromaser, and thus improve the chances for realizing 
connections between the detection statistics and the field variables in a working 
micromaser. Possibilities include two-atom events, spontaneous emission, veloc- 
ity averaging, and suppression of perfect inversion by thermal photons, among 
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others. Also important is knowing which of these effects is the most significant. 
We plan to pursue this direction of research in the future. 
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Figure 1: Plots (a) and (b) show the atomic inversion for A^ox — 7 and A^ox = 10 
respectively, for comparison to Figs. 3{a) and 3(/3) of Ref. [|lj. Here i' = 0.054 
and r] = 40%. The vertical dotted lines indicate the same trapping states as in 
Fig. 3 of Ref. jl]. The dashed lines show the effect of a 3 /xs uncertainty in tint 
with gaussian distribution. 
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Figure 2: Steady state Fano-Mandel function Qa(oo) for detections of excited 
atoms for iVex = 7 and v = 0.054. The result is independent of the detector 
efBciencies when scaled by l/r?A- 




Figure 3: Steady-state Fano-Mandel function (5b(oo) for detections of de-excited 
atoms for A'ex = 7 and u = 0.054. The result is independent of the detector 
efficiencies when scaled by I/t/b- 
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iint (fis) 

Figure 4: Fano-Mandel function QB{t) for detections of de-excited atoms for 
comparison to Fig. 4 of Ref. [Q. Here A'ox — 7 and v = 0.054. The long-dashed 
line shows Quit = 4/7), the dashed Hue shows Quit = I/7), and the solid line 
shows the result of averaging Qb {t) for times between one and four cavity decay 
times. The vertical dotted lines indicate the same trapping states as in Fig. 4 
of Ref. §. 
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Figure 5: Mean number {ua) of successive detections of excited atoms for Nex = 
7 and v = 0.054. The curves are for ij = 100% (solid hne), ij = 40% (long dashed 
line), and 77 = 10% (dashed line). The vertical dotted lines indicate the trapping 
states for no = 0, ... , 4. 
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Figure 6: Mean number (ne) of successive detections of de-excited atoms for 
iVex = 7 and v = 0.054. The curves are for r] = 100% (solid line), r] = 40% (long 
dashed line), and -q = 10% (dashed line). The vertical dotted lines indicate the 
trapping states for no = 0, . . . , 4. 




Figure 7: Mean number of successive detections (of any type) normalized to the 
uncorrelated value for iVex = 7 and v = 0.054. The curves are for r] = 100% 
(solid line), r] = 40% (long dashed line), and r] = 10% (dashed line). 
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Figure 8: Mean waiting time (^a^b) until the next detection of a de-excited 
atom after an initial detection of an excited atom (scaled by the rvye) for A'ox = ^ 
and V = 0.054. The curves are for -q = 100% (solid line), r] = 40% (long dashed 
line), and r] = 10% (dashed line). The vertical dotted lines indicate the trapping 
states for no = 0, ... , 4. 




Figure 9: Moan waiting time {Ib^a) until the next detection of an excited 
atom after an initial detection of a de-excited atom (scaled by ri]A) for A^ex = 7 
and u = 0.054. The curves are for r] = 100% (solid line), r] = 40% (long dashed 
line), and r] = 10% (dashed line). The vertical dotted lines indicate the trapping 
states for no = 0, ... , 4. 
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Figure 10: Mean waiting time (^a^b) until the next detection of a de-excited 
atom after an initial detection of an excited atom normalized to the uncorrelated 
value for TVcx = 7 and v = 0.054. The curves are for t] = 100% (solid line), 
r] = 40% (long dashed hne), and r/ = 10% (dashed line). 




Figure 11: Mean waiting time {I-q^a) until the next detection of an excited atom 
after an initial detection of a dc-cxcited atom normalized to the uncorrelated 
value for A^ex = 7 and u = 0.054. The curves are for r] = 100% (solid line), 
T} = 40% (long dashed line), and r} = 10% (dashed line). 
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